#################################################################################
#############################Market Potential####################################

pm<-function(data,W,mass.var=NULL){
stopifnot(nrow(W)==ncol(W))
if(mass.var=='NULL'){
	stopifnot(is.vector(x)==TRUE)
	x<-data
	}else{
	x<-data[,mass.var]
	}
 pm<-x/W
  for(i in 1:nrow(W)){
	pm[pm[,i]==Inf,i]<-0
  	pm[is.na(pm[,i]),i]<-0
	}

 pm<-colSums(pm)
 pmR<-(pm/max(pm))*100
	z<-list('pmA'=pm,'pmR'=pmR)
	class(z)<-'pm'
	return(z)
 }


getwd()

#save(pm,file='PM.q')



